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NOTATIONS 


a  Radius  of  dust  particles,  centimeters 

z  Number  of  electronic  charges  on  dust,  a  positive  or  negative 

integer 

e  Electronic  charge 

T  Absolute  temperature 

"C  Mean  thermal  speed  of  ions  or  electrons. 

For  ions,  assume  a  mass  of  50  amu. 

"*  1 

Therefore  C  +  -  2.6  x  10  cm  s  .  For  electrons,  *  300  C  + 

F  Geometric  modification  shape  factor 

Ng  Number  density  of  electrons,  centimeter" 3 

N+  Number  density  of  positive  ions,  centimeter"3 

N_  Number  density  of  negative  ions,  centimeter" 3 

Np  Number  density  of  dust,  centimeter" 3 

vz  Frequency  of  capture  of  positive  ions  by  dust  particles 

with  a  charge  z,  second'1 

V®  Frequency  of  capture  of  electrons  by  dust  particles  with 

a  charge  z,  second'1 

v”  Frequency  of  capture  of  negative  ions  by  dust  particles 

'  with  a  charge  z,  second"1 

vph  Photodetachment  frequency,  second"1 

Pz  Fraction  of  dust  particles  which  carry  a  charge  z 

0 

kp  Coefficient  of  attachment  of  an  electron  to  a  single 

dust  particle,  centimeter3  second"1 

kp  Coefficient  of  attachment  of  a  positive  ion  to  a  single 

1  dust  particle,  centimeter3  second"1 

K®  Effective  attachment  coefficient  of  electrons  in  a 

u  distribution  of  dust  particles,  centimeter3  second"1 

Kp  Effective  attachment  coefficient  of  positive  ions  in 

a  distribution  of  particles,  centimeter3  second"1 
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1.  INTRODUCTION 


The  D  region  is  the  lowest  portion  of  the  ionosphere  and  nominally  encompasses 
the  altitude  range  of  60  to  90  km.  It  represents  a  transition  region  for  the 
ionospheric  plasma  composition.  Above  80  to  85  km,  daytime  composition  is 

limited  to  molecular  ions  such  as  N0+  and  0^;  while  below  80  to  85  km,  the 

simple  molecular  ions  enter  into  clustering  reactions  which  convert  the 

positive  ion  composition  into  one  dominated  by  water  clusters.  Electrons  also 

undergo  attachment  processes  in  the  D  region  to  form  negative  ions;  by  day 

this  attachment  appears  to  take  place  below  70  to  75  km.  By  night  it  occurs 

as  high  as  90  km,  but  there  are  still  many  unanswered  questions  concerning  the 

negative  ion  chemistry. 


The  sources  of  ionization  within  the  D  region  are  the  "tail  end"  of  the  solar 

x-ray  spectrum,  mainly  at  wavelengths  less  than  1  nm  and  in  the  3-  to  10-nm 

range,  and  two  strong  solar  emission  lines,  Lyman  beta  at  102.6  nm  which 

ionizes  02  and  Lyman  alpha  at  121.6  nm  which  ionizes  NO.  Ionization  from 

these  sources  occurs  principally  above  75  km.  Below  70  km  galactic  cosmic 
'  rays  are  the  main  (though  small)  source  of  ionization.  At  higher  latitudes 

charged  particle  precipitation  and  associated  bremsstrahlung  x-rays  may 
occasionally  be  the  dominant  ionization  source. 

The  deposition  of  micrometeors,  at  the  altitude  of  the  lower  portion  of  the  E 
region  and  upper  portion  of  the  0  region,  produces  a  concentration  of  metallic 
ion  layers  at  90  to  110  km  and  occasional  occurrences  of  sporadic  E 
reflection.  Estimates  have  shown  that  this  micrometeor  flux  may  be 
sufficiently  large  to  serve  as  source  material  for  the  mesospheric  scattering 
layers  and  nocti lucent  clouds  which  are  often  observed  in  the  80  to  85  km 
range.  1 

Certain  D  region  phenomena  have  been  observed  which  have  no  ready  explanation 
in  terms  of  conventional  ion  chemistry;  for  example,  electron  density  profiles 
have  occasionally  shown  layers  of  sharply  decreased  concentration  in  the  30  to 
90  km  range,  which  could  be  interpreted  as  due  to  particulate  layers  which 
have  the  ability  to  attach  electrons.  2  This  ability  would  suggest,  therefore, 
that  there  is  a  large  group  of  primarily  negatively  charged  species,  either 


‘0.  M.  Hunten,  R.  P.  Turco,  and  0.  B.  Toon,  1980,  "Smoke  and  Dust  Particles  of 
Meteoric  Origin  in  the  Mesosphere  and  Stratosphere,"  J  Atmos  Sci,  37:1342. 

*A.  Pedersen,  J.  Troin,  and  J.  A.  Kane,  1970,  "Rocket  Measurements  showing  a 
Removal  of  Electrons  Above  the  Mesopause  in  Summer  at  High  Latitude,"  Planet 
Space  Sci ,  18:945-947. 
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ions  or  particulates,  present  in  the  upper  D  region.  3  '  And  indeed, 

conductivity  measurements  throughout  the  stratosphere  and  mesosphere  give 
evidence  of  light  and  "massive"  ions,  which  may  be  interpreted  to  be  changed 
particulates.  D-region  electron  densities  observed  during  some  solar  eclipses 
have  indicated  a  surprisingly  rapid  electron  density  decrease  which  cannot  be 
explained  on  the  basis  of  current  knowledge  of  electron-ion  nates  of 
recombination.  5  3 

With  these  examples  of  anomalous  behavior  in  mind,  the  purpose  of  t^'s  rP,,'cr‘ 
is  to  establish  properties  of  particulates  necessary  for  tso  c' c- 

predictive  computer  codes.  These  properties  are  electron-parti cu1 at-  ann 
particulate  collision  frequencies,  attachment  and  reconvvnati-jn  rat  -, 
electron  photodetachment  rates.  To  achieve  this  result,  the  cozee'C:‘?''“ ' 
be  examined  on  three  levels:  (1)  Coefficients  will  be  der-'vcd  for  s ~ •  ; , •  ■ 
particles  impacting  with  electrons  or  ions.  (2)  An  "effective  cr  >. ’ 

attachment  coefficient  shall  be  derived  which  pertains  to  a  population  of  -t:  : 
mixed  with  ambient  electrons  and  positive  ions.  This  derivation  fo1! ows  t 
work  of  Parthasarathy  and  Ra i  6  and  Parthasarathy 7  but  has  certain  Iviit-*-- 
which  shall  be  examined.  This  treatment  is  instructive  since  it  derives  t"> 
charge  loading  on  a  population  of  dust  and  outlines  the  sensitivity  of  the 
attachment  coefficients  to  various  parameters.  (3)  Particulates  w' '  '  n 
included  and  evaluated  as  electron  or  ion  attachment  sites  within  a  ful 1  time- 
dependent  chemical  code  which  includes  all  other  competitive  processes  c 4 
region  chemistry. 


3M .  G.  Heaps,  F.  E.  Niles,  and  R.  D.  Sears,  1978,  Modeling  the  Ion  Chemistry 
of  the  D-region:  A  Case  Study  Based  Upon  the  19613  TotaT  Solar  Tcllpse, 
ASl-TR-0015,  US  Army  "Atmospheric  Sciences  Laboratory,  White  Sands  MissTTe 
Range,  NM. 

-M.  G.  Heaps  and  J.  M.  Heimerl ,  1980,  "The  Quiet  Midlatitude  D-region:  A 
Comparison  Between  Modeling  Efforts  and  Experimental  Efforts,"  J  Atmos  Terr 
Phys,  42:733. 

5E.  A.  Mechtly,  C.  F.  Sechrist,  Jr.,  and  L.  G.  Smith,  1972,  "Electron  Loss 
Coefficients  for  the  0-region  and  the  Ionosphere  from  Rocket  Measurements 
During  the  Eclipses  of  March  1970  and  November  1966,"  J  Atmos  Terrs  Phys, 
34:641 

‘R.  Parthasarathy  and  D.  B.  Rai,  1966,  "Effect  of  Meteoric  Oust  on  the 
Effective  Recombination  Coefficient  in  the  Lower  Ionosphere,"  Radio  Sci, 
1:1401. 

7R.  Parthasarathy,  1976,  "Mesopause  Oust  as  a  Sink  for  Ionization,"  J  Geophys 
Res,  81:2392. 
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2.  PROPERTIES  OF  MESOSPHERIC  PARTICLES 


Noctilucent  clouds  (NIC)  have  long  provided  the  most  direct  and  visual 
evidence  that  particulates,  most  likely  pure  ice  or  ice  coated  dust,  tend  to 
layer  at  the  summer  mesopause.  The  rare  observation  of  NLC  during 
stratospheric  warmings  indicates  layering  in  winter  also.8  0G0-6  satellite 
observations’  indicate  that  particulate  layers  near  the  82  km  upper  bound  of 
the  mesosphere  may  be  a  persistent  high  latitude  feature,  while  global  lidar 
observations  show  a  sodium  layer  in  this  same  region.  Witt10  and  Fogle  and 
Rees11  have  calculated  NLC  particulate  sizes  to  be  of  the  order  of  O.lum  to 
0.2um  in  radius  with  number  densities  of  the  order  of  1  cm"’.  Donahue  et  al * 
also  estimate  that  0.13un  radius  particles  with  number  densities  of  15  to  40 
cm"3  would  successfully  explain  their  scattering  observations.  An  earlier 
work  by  Divari12  based  upon  rocket  samples  states  that  particulate 
concentrations  in  the  mesosphere  are  an  inverse  function  of  size  down  to 
0.05um.  Also  model  calculations  by  Hunten  et  al  1  indicate  a  similar  inverse 
size  and  concentration  relationship  in  the  80-  to  90-km  altitude  range  for 
particles  less  than  O.OOlum. 

There  are  three  alternatives  for  the  origin  of  NLC  particulates:  (1)  They  are 
terrestrial  and  ascend  by  the  action  of  turbulence,  (2)  they  are  formed  from 
in-situ  ion  clusters,  or  (3)  they  are  meteoritic. 

The  first  possibility  seems  unlikely  considering  the  insufficient  magnitude  of 
the  eddy  diffusion  coefficient  in  the  stratosphere  and  mesosphere  (see,  for 
example,  Hunten13).  As  for  the  second  possibility,  Witt14  proposed  that  water 


9R.  A.  Hamilton,  1964,  "Observations  of  Noctilucent  Clouds  Near  Midwinter," 
Met  Mag  Lond,  93:201. 

’T.  M.  Donahue,  B.  Guenther,  and  J.  E.  Blamont,  1972,  "Noctilucent  Clouds  in 
the  Daytime:  Circumpolar  Particulate  Layers  Near  the  Summer  Mesopause," 
Atmos  Sci ,  29(6) : 1 205 . 

I4G.  Witt,  1960,  "Polarization  of  Light  From  Noctilucent  Clouds,"  J  Geophys 
Res,  65:925-933. 

‘‘B.  Fogle  and  M.  Rees,  1972,  "Spectral  Measurements  of  Noctilucent  Clouds," 
Geophys  Res.  77:720. 

12G.  Divari,  1964,  "The  Dust  Concentration  in  the  Upper  Layers  of  the  Earth’s 
Atmosphere,"  Geomagnetism  and  Aeronomy,  4:688-692. 

*D.  M.  Hunten,  R.  P.  Turco,  and  0.  B.  Toon,  1980,  "Smoke  and  Dust  Particles  of 
Meteoric  Origin  in  the  Mesosphere  and  Stratosphere,"  J  Atmos  Sci,  37:1342. 

13D.  M.  Hunten,  1975,  Vertical  Transport  in  Atmospheres,  Atmospheres  of  Earth 
and  the  Planets,  D.  Reidel  Publishing  Co.,  Oordrecht-Hoi land,  pp  5 1-/2. 

14G.  Witt,  1969,  "The  Nature  of  Noctilucent  Clouds,"  Space  Res,  IX :  157-169 , 
North  Holland  Publishing  Co. 
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cluster  ions  in  the  mesosphere1’  16  serve  as  condensation  or  sublimation 
nuclei.  Such  a  clustering  mechanism  has  been  investigated  for  aerosols  in 

general.17  19  However,  the  mechanism  whereby  clusters  of  ten  or  fewer 
molecules  could  grow  or  combine  to  form  the  size  of  particulates  nominally 
present  in  the  D  region  (hundreds  or  thousands  of  molecules)  has  not  yet  been 
identified.  Assisted  by  a  model  calculation,  Reid19  has  argued  that  if  a 

substratum  of  small  particulates  of  nominal  O.OOlum  size  exists,  then  the 
growth  of  ice  needles  of  a  detectable  size  (0. 13um)  is  possible  under 

constraints  of  mesospheric  temperatures  and  low  water  vapor.  The  idea  tnat 
meteoritic  dust  particles  could  be  sufficiently  numerous  to  form  layers  dense 
enough  for  NLC  also  has  been  theoretically  supported  by  Chapman  and 
Kendall.  20  Hunten  et  al  1  showed  that  a  nominal  meteor  flux  of  lQ"lf  o  c"'": 
s'1  is  sufficient  to  provide  a  source  of  particulates  with  sizes  of  the  cwr 
of  0.001am  and  number  densities  of  10!  cm"'.  Also  the  source  of  particles  may 
be  multiple  and  seasonally  dependent  as  suggested  by  data  of  Clemesha  ot 

al.21  Thus,  mesospheric  particulate  sizes  and  their  number  densities  can 
apparently  range  from  an  upper  bound  of  0.1am  radius  with  number  densities  or’ 
the  order  of  10  cm"  3  or  less,  to  a  lower  bound  of  O.OOlum  sized  particulates 
with  number  densities  on  the  order  of  103  cm'3.  This  range  of  concentrations 
and  sizes  will  be  used  as  representative  in  the  subsequent  calculations  and 
modeling  in  this  report. 


15R.  0.  Narcisi  and  A.  D.  Bailey,  1%5,  "Mass  Spectrometri c  Measurements  of 
Positive  Ions  at  Altitudes  From  64  to  112  Kilometers,1  J  Geophys  Res,  ’70:.?0F  • 

1 6F .  Arnold  and  W.  Joos,  1979  ,  "Rapid  Growth  of  Atmospheric  Cluster  Ions  at 
the  Cold  Mesopause,"  Geophys  Res  letters,  6:763-766. 

17A.  W.  Castleman,  Jr.,  1974,  "Nucleation  Processes  and  Aerosol  Chemistry," 
Space  Sci  Rev,  15:547- 589 . 

l,A.  W.  Castleman,  Jr.,  1979,  "Nucleation  and  Molecular  Clustering  About 
Ions,"  Adv  in  Colloid  and  Interface  Sci,  10:73-128. 

1 ’G.  C.  Reid,  1975,  "Ice  Clouds  at  the  Summer  Polar  Mesopause."  J  Atmos  Sci. 
32:523. 

J0S.  Chapman  and  P.  C.  Kendall,  1965,  "Noctilucent  Clouds  and  Thermospheric 
Dust:  Their  Diffusion  and  Height  Distribution,"  Quart  J  Roy  Met  Soc,  91:115. 

‘D.  M.  Hunten,  R.  P.  Turco,  and  0.  B.  Toon,  1980,  "Smoke  and  Dust  Particles  of 
Meteoric  Origin  in  the  Mesosphere  and  Stratosphere,"  J  Atmos  Sci,  37:1342. 

21B.  R.  Clemesha,  V.  W.  J.  H.  Kirchoff,  D.  M.  Simonich,  and  H.  Takahashi, 
1978,  "Evidence  of  an  Extra-Terrestrial  Source  for  the  Mesospheric  Sodium 
Layer,"  Geophys  Res  Letters,  5:873-877. 
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3.  ELECTRON  AND  ION  ATTACHMENT  TO  SINGLE  DUST  PARTICLES 

The  first  quantitative  steps  towards  examining  the  properties  of  dust  in  the 
presence  of  charged  particles  were  taken  in  1960  by  Natanson,22  who  formulated 
a  kinetic  approach  to  the  charging  of  dust  by  electrons  and  ions. 
Parthasarathy  and  Rai  6  used  Natanson's  approach  to  derive  attachment 
coefficients  for  dust  particles  with  electrons  and  positive  ions  in  the 
mesosphere.  This  work  was  extended  by  Parthasarathy 7  to  include  a 

probabilistic  formulation  of  the  basic  equilibrium  equations  including  an 
explicit  photodetachment  term. 

A  dust  particle  is  charged  when  an  electron  or  ' on  collides  with  and  is 
captured  by  the  particle.  The  capture  r'requency,  v,  equa7  to  the  collision 
frequency  (with  a  collision  efficiency  of  1),  may  be  expressed  as 


.  =  ta  N  i  F  is  )  , 


where  ^a2  as  the  geometric  c-'oss  section  of  the  particle,  N  is  electron  or  ion 
density,  C  is  mean  velocity,  and  F  is  a  capture  rate  modification  function  or 
"shape  factor'  7  which  modifies  the  geometric  cross-section  and  whose  value 
depends  upon  the  resident  charge  of  the  particle  and  its  radius.  This  shape 
factor,  F,  expresses  the  attractive  (or  repu1sive)  effects  of  charged  dust 
particles  to  approaching  electrons  or  ions  (figure  1). 

Several  points  should  be  noted.  First,  the  F  factor  for  uncharged  dust 
particles  enhances  the  geometric  cross  section,  an  effect  which  increases  as 
the  particle  shrinks.  This  effect  is  due  to  the  setting  up  of  an  image  charge 
within  the  particle  as  the  electron  or  ion  approaches  its  surface.  Second, 
powerful  repulsive  effects  lead  to  a  practical  limit  of  one  charge  for 
particles  of  approximately  O.Olum  in  size  or  less.  For  such  small  particles, 
recombination  or  neutral ization  becomes  greatly  enhanced  as  the  particle 
becomes  charged.  Only  for  particles  larger  than  O.lum  does  multiple  charging 
become  probable. 


22G.  L.  Natanson,  I960,  "On  the  Theory  of  the  Charging  of  Amicroscopic  Aerosol 
Particles  as  a  Result  of  the  Capture  of  Gas  Ions,"  Sov  Phys  Tech  Phys, 
(English  Translation),  5:538. 

6R.  Parthasarathy  and  D.  B.  Rai,  1966,  "Effect  of  Meteoric  Dust  on  the 
Effective  Recombination  Coefficient  in  the  Lower  Ionosphere,"  Radio  Sci , 
1:1401.  . 

’R.  Parthasarathy,  1976,  "Mesopause  Dust  as  a  Sink  for  Ionization,"  J  Geophys 
Res,  81:2392 
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CAPTURE  RATE  MODIFICATION  FUNCTION,  (F) 


0  12  3  4 

ABSOLUTE  CHARGE  Z 


A  plot  of  the  capture  rate  modification  function  (F)  against  the 
absolute  charge,  Z,  on  a  dust  particle  of  radius  O.lum  to 
O.OOlum.  The  value  of  the  capture  rate  modification  function 
expresses  the  ratio  of  the  "elective"  capture  cross  section  of  a 
dust  particle  to  the  simple  geometric  cross  section  and  is 
dependent  upon  electrostatic  effects. 


12 


The  dust  charge  attachment  coefficient,  kD,  for  an  isolated  particle  is 
defined  as 


kp  =  v/N  =  ua  C  F  (cm  s’  )  , 


(2) 


where  k0  is  a  function  of  particle  size,  resident  charge,  electron  or  ion 
mass,  and  velocity,  kp  is  to  be  treated  as  a  standard  chemical  rate  constant, 
where  the  term  in  the  set  of  continuity  equations  becomes 
kp  [Np]  [Ng]  or  kp  [Np]  [N+  ]  .  Table  1  lists  representative  values  of  the 
coefficients  kp  and  kp  for  a  nominal  mesospheric  temperature  T  =  15CTK. 

Points  which  are  apparent  from  table  1  are:  (1)  the  smaller  the  dust,  the 
more  rapidly  the  electron  attachment  coefficient  decreases  with  increasing 
negative  charge,  (2)  the  larger  the  particle,  the  greater  tendency  it  will 
have  to  attach  electrons  and,  hence,  provide  an  electron  sink.  A  single  large 
particle  (that  is,  radius  >0.1um)  is  therefore  as  effective  as  several  smaller 
particles  in  perturbing  the  charge  balance. 


TABLE  1.  kp  A HD  kp ,  ELECTRON  AND  POSITIVE  ION  ATTACHMENT 

COEFFICIENTS  FOR  A  SINGLE  OUST  PARTICLE  AT  T  =  150°K 


Radius  (urn) 

Charge  (e) 

kp  (cm5  s' ’) 

kp  (cm3  s'1) 

Z 

0.001 

+1 

2 . 6  ( -  5) 

1 . 3  ( -  3  3 ) 

0 

3 . 3  ( -  6 ) 

1.0(-8) 

-1 

4 . 2  ( -  3 1 ) 

8 . 4 ( -8) 

-2 

2 . 3 ( -37) 

1 . 6  ( -  7 ) 

0.01 

+  1 

2. 8( -4) 

7 . 3( -1 1 ) 

0 

1 . 2( -4 ) 

3.9(-7) 

-1 

2 . 3  ( -  7 } 

9.1(-7) 

-2 

l.K-10) 

1 .  7  ( -6 ) 

0.1 

+  1 

5. Q( -3) 

l.K-5) 

0 

5 . 5  ( -  3 ) 

1 . 7( -5) 

-1 

3 . 4  ( -  3 ) 

1 . 5( -5) 

-2 

1 . 3  (  -  3 ) 

2 . 4  ( -  5 ) 

-3 

5 . 2  ( -  4 ) 

3 . 2  ( —  5) 

-4 

2.1(-4) 

4 . 0  ( -  5 ) 

-5 

8.3(-5) 

4.9(-5) 

-6 

3. 3( -5) 

5 . 7  { — 5 ) 

-7 

1 . 3 ( -5) 

6 . 5  ( -  6 ) 

-8 

5 . 2  ( -  6 ) 

7.4( -6) 
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A  further  insight  may  be  provided  by  a  plot  (finu>*e  ?)  comparing  the  magnif  y 

of  the  attachment  coefficients,  k(, ,  for  z  .  0,  and  +1  charged  dust,  as 

function  of  particle  radius,  and  the  elect.ron-ion  recombination  coef t  i l i •■■■t 
for  N0  +  and  water  cluster  ions.  This  figure  shows  that  although  normally  th. 
most  important  loss  mechanism  for  electrons  is  recombi nati on  to  cluster  ions, 
electron  attachment  to  particulates  couM  be  dominant  if  the  dust  ny» 
density  were  of  the  same  magnitude  as  the  ion  number  density  and  'a'  -;ust  «  ■■ 
positively  charged,  (b)  dust  were  neutral  but  '  anger  thin  vt  n,  o*  '  *’• 

dust  particle  had  a  single  negative  change  ,v  :  ex"eeded  .  : 

The  inferences  of  dust  size  distribution  *h<.  mesesph^'-o  oi^  fa*  ♦  ►  • 
could  be  much  dust  at  sizes  of  J.Cl  -.or  or  smal lor.  In  f'at  „  *■  ■ 

probability  of  mesospheric  dust  affect  He  the  "'-''•eqioi  rr,v  ge  ra  ■'-■■■  >  , 
be  considerable  through  the  formation  oe  an  •  *•-<  n  s'-k. 

4.  EFFECT  I VE  CHARGE  ATTACHMENT  COEFFICIENTS 

4.1  Derivation 

The  previous  section  explicitly  derived  the  charge  attachment  frequency  ‘  •  . 

isolated  dust  particle  of  a  given  size  and  a  known  resident  charge.  --„r.v 
in  the  D  region,  a  population  of  dust  parti  les  will  tend  to  puss*  -  ■ 

distribution  of  charge  which  will  be  estab1  '  shod  by  the  competitive  at  fa  '■  • 

and  detachment  of  electrons  and  ions.  Darthasarathy  and  Pa- 
Parthasarathy 7  initially  addressed  the  calculation  of  the  attuf*  ■ 
coefficients  for  dust  "en  masse."  In  this  section  we  analyze  the  limitaf  ■ 
of  this  approach  but  show  that  <t  provides  knowledge  of  the  resident  v  :• 
distribution  function  of  dust  for  various  anh'oni  conditions.  c hotoneta'  n- •• 
is  also  included  in  this  formulation  a  a->  explicit  term,  and  su  i  *  ib‘  • 

experimental  values  are  examined. 

Parthasarathy 7  stipulates  that  the  equilibrium  condition  of  dust  mixed  with 
charged  species  is  considered  to  be  a  sf* * st ’ tally  stationary  di stri but •  m . 

This  condition  is  formulated  by  stating  that  the  frequency  at  which  the 

fraction  of  dust  particles  with  a  charge  z  is  being  depleted  by  ion  capture, 
electron  capture,  and  photodetachment  is  equal  to  that  frequency  at  which  the 
fraction  of  dust  particles  with  z  +  1  charges  are  acquiring  electrons  or 
negative  ions  plus  the  frequency  at  which  the  fraction  of  dust  particles  with 
z  -  1  charges  are  acquiring  positive  charges  and  suffering  photodetachment  cf 
electrons  (after  Parthasarathy 7). 


‘R.  Parthasarathy  and  D.  B.  Rai,  1966,  "Effect  of  Meteoric  Dust  on  the 
Effective  Recombination  Coefficient  in  the  Lower  Ionosphere,"  Radio  Sci, 
1:1401. 

7R.  Parthasarathy,  1976,  "Mesopause  Dust  as  a  Sink  for  Ionization,"  J  Geophys 
Res,  81:2392. 
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fcp, DUST-ELECTRON  ATTACHMENT  COEFFICIENT  (cm 


ELECTRON-ION  RECOMBINATION  COEFFICIENT  ( 


Expressed  algebraically. 


Pz  (v+z  + 


ve  + 
Z 


vph^  Pz+1 


:Vi  * 


z  +  1 


+  P 


z-1 


Vi  ♦ 


(3) 


By  setting  Pz  =  0  for  a  large  value  of  charge  (say  z  =  -9,  meaning  that  no 
dust  has  9  or  more  electrons  on  it),  and  solving  the  other  Pz's  in  terms  of 
P_g  and  so  forth,  a  recursion  relation  may  be  obtained.  Then,  after 
normalizing  with  iP  -  1,  solve  for  each  P  Neglecting  negative  ions  for  t^e 

daytime  0  region,  the  effective  attachment  coefficients  are  then  def-'ned  as: 


=  i  rP  v+  =  i  rP  (ve  -  v  J 
D  n  ,  z  z  N  i  z  ph 
+  e 


v+  -  1  rP  e  l?z  •  V’Ph 
D  N  z  z  N 
e  e 


The  attachment  coefficient  for  electrons,  Kr^,  is 


,e  ePz”2 
'D  = 


Substituting  equation  (6)  into  (5)  and  remembering  sP2  -  i. 


K 


0 


(7) 


that  is,  the  effective  attachment  coefficient  for  positive  ions  is  equal  to 
the  effective  attachment  coefficient  for  electrons  diminished  by  the 
photodetachment  frequency  divided  by  the  ambient  electron  density. 

This  formulation  incorporates  the  following  assumptions  and  restrictions  which 
limit  i cs  applicability  to  D  region  time  dependent  chemical  codes: 

a.  This  formulation  is  explicitly  stationary.  An  equilibrium  is  assumed 
between  the  rate  of  attachment  and  detachment  of  ions  and  electrons  for  the 
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derivation  of  the  P2  fractions  and  the  coefficients.  As  such,  these 

coefficients  are  unsuitable  for  use  in  time  dependent  programs  designed  for 
rapidly  varying  atmospheric  conditions  such  as  eclipse  times,  sunrise  and 
sunset,  or  nuclear  burst  phenomena. 

b.  It  is  assumed  that  Ne  =  rN+.  This  assumption  of  charge  neutrality 

among  electrons  and  ions  is  equivalent  to  the  assumption  that  the  number  of 
particulates  is  too  small  to  perturb  the  charge  balance,  which  may  be  an 
unrealistic  assumption  for  the  mesosphere.  Therefore,  since  the  derived  P2 

values  are  func;  ins  of  the  Ne  or  rN+  used,  the  calculation  of  new  °2 ' s  and 

the  effective  coefficients  would  have  to  be  handled  iteratively  in  a  computer 
code  each  time  new  Np  and  cN+  were  computed. 

c.  This  formulation  (equation  3)  neglects  competitive  chemical  attachment 
processes  such  as  the  recombination  reactions  between  positive  ions  and 

0  + 

electrons  which  would  change  °2  values  and  Krj  and  . 

0  + 

Within  these  restrictions  the  values  of  KQ  and  derived  by 

Parthasarathy 7  are  useful  for  steady  state  chemical  codes  with  a  comparatively 

small  number  of  particulates.  But  if  no  restriction  is  to  be  placed  upon  the 
size  of  the  dust  population  or  upon  the  concentrations  of  electrons  or 
positive  ions,  then  a  set  of  stepwise  reactions  utilizing  the  coefficients  for 
individual  dust  in  various  charge  states  is  more  appropriate  (section  2). 
However,  a  benefit  of  the  approach  by  Parthasarathy  is  to  predict  the  range 
and  stability  of  resident  charge  loading  expected  on  a  population  of  dust 

under  various  ambient  conditions.  Thus,  the  number  of  piecewise  reactions 

sufficient  for  dust  of  a  given  size  may  be  predicted  for  use  in  the  chemical 
reaction  code. 

4.2  Choice  of  Photodetachment  Coefficient  for  Dust 

The  photodetachment  rate  for  the  dust  particle,  vph,  is  a  key  and  governing 

parameter  in  these  computations.  Photodetachment  yields  of  many  materials  are 
difficult  to  measure  in  a  laboratory,  and,  in  turn,  whether  such  measurements 
would  be  valid  for  mesospheric  conditions  is  difficult  to  assess.  The  best 
that  may  be  done  is  to  assume  a  likely  composition  fo>-  dust  and  establish  a 
plausible  range  for  the  photoelectric  yield  over  the  incident  solar 
wavelengths  which  exist  in  the  mesosphere. 

Photoelectric  yield  of  a  given  material  is  a  function  of  many  parameters  such 
as  surface  roughness,  angle  of  incidence  of  photons,  temperature,  and  possible 
adsorption  of  gases  on  the  surface.  For  example,  photoelectric  yield  will 


’R.  Parthasarathy,  1976,  "Mesopause  Dust  as  a  Sink  for  Ionization,"  J  Geophys 
Res,  81:2392. 


increase  with  greater  roughness  and  angle  of  incidence.2’  Adsorption  of  gases 
such  as  molecular  oxygen  or  hydrogen  onto  the  surface  may  reduce  the  yield  and 
shift  the  cutoff  wavelength  further  into  the  ul traviolet;  2 '  whereas  very  small 
amounts  of  gases  may  increase  yield  and  shift  the  cutoff  wavelength  to  longer 
wavelengths  through  the  formation  of  surface  "double  layers"  which  a^low  the 

creation  of  surface  states  about  0.4  ev  below  the  Fermi  level.2* 

If  meteoritic  material  from  galactic  space,  possibly  ice  coated,  is  assumed  t- 
be  the  constituent  of  the  dust  particles,  then  the  photodetachment.  values  fr- 
moon  dust,25  nickel,23  or  ice6  would  be  reasonable  options,  "resunab!  v . 
region  meteoritic  depositions  would  bc>  of  similar  material  to  1  ->r  c-j-* 
dust.  Will’s  et  al25  determined  the  photoelectric  yield  of  sarv'es  or  ’c  > 
dust  from  the  Apollo  14  and  15  mi ssi ons  "•  a<  ipri-'-;eations  r v  Tn 

were  made  at  10"’  to  10”’  torr,  a  more  abso’ute  vacuum  V  an  in  the  ;*  cm  y 
five  orders  of  magnitude,  impurities  ir  the  s-u-pl es  we^e  baked  out  for  ’ 
at  150°C .  Auger  electron  spectroscopy  ?ng  hem’cal  jna1  /«.’  s  j'?'  dr  ,J 
composition  of:  48  percent  of  S’Oo,  IF  percen»  of  11 

C  c 

10  percent  of  FeO,  etc.,  which  agrees  with  other  determinations  '.nuop  ■■ 
for  si  1 icate-type  meteors.  26  Peak  emission  of  10" :  electrons  .>•?•*  :  hf>i ■■  ■: 
occurs  at  800  angstroms,  with  smaller  va1  ups  of  in"6  at  >!«•'.  .vastrue.-  >.  ■ - 
10* 1 2  at  600  angstroms.  ror  solar  flux  at  the  top  cl  th'-'  earth’s  atmospr'-*  •• 
the  composite  yield  is  2.8  x  10’  ( wa 2 )  photoel  ectrens  per  second  for  l;jn.;-- 
surface  material,  where  ita 2  represents  the  surface  area  in  square  centi¬ 
meters.  Attenuation  of  solar  flux  in  the  Schumann-Runge  Rands  and  Lyina*'  i  '.r' 
the  mesosphere  reduces  this  photodetachment  frequency  of  neutral  duct  at  80  km. 
to  about  1  x  10’  ( wa 2 1  s’1.  Since  iron  and  nickel  meteorites  have  been 
commonly  reported,  26  the  photoelectric  yield  for  elemental  nickel  was  obtained 
in  Weissler.23  When  the  appropriate  solar  flux  is  used  for  Lyman  t  and  the 

Schumann-Runge  8ands  attenuation  at  80  km  is  accounted  for,  the  resultant 

photodetachment  for  nickel  is  1.5  x  10"  (?a2'  s'1.  Values  of  photodetachment 
for  ice6  are  1.9  x  10 *  ( *a 2)  s"‘.  To  bracket  these  previous  determinations 
for  different  materials  and  to  allow  for  the  considerable  uncertainty  which 
still  remains,  the  following  values  of  photodetachment  frequency  are  chosen  to 


2IG.  L.  Weissler,  1956,  "Photoionization  in  Gases  and  Photoelectric  Emission 
from  Solids,"  in  Handbuch  der  Physik,  Vol  XXI,  Springer-Verl ag,  304-33?. 

2*B.  Feuerbacher  and  B.  Fitton,  1972,  "Photoemission  from  Surface  States  on 
Tungsten,"  PhysRey  Letters,  29(121:786. 


25R.  F.  Willis,  M.  Anderegg,  B.  Feuerbacher,  and  B.  Fitton,  "Photoemission  and 
Secondary  Electron  Emission  from  Lunar  Surface  Material,"  in  Photon  and 
Particle  Interactions  with  Surfaces  in  Space,  Proceedings  of  the'  5th  Eslab' 
Symposium,  Noordwijk,  the  Netherlands,  26-29  September  1972.  R.  1.  L.  Grand, 
Ed.,  D.  Reidel  Publishing  Co.,  Dordrecht-Hol 1  and/Boston  -  USA. 

‘R.  Parthasarathy  and  0.  B.  Rai,  1966,  "Effect  of  Meteoric  Oust  on  the 
Effective  Recombination  Coefficient  in  the  Lower  Ionosphere,"  Radio  Sci. 
1:1401.  - — 

2,J.  A.  Jacobs,  R.  D.  Russell,  and  J.  T.  Wilson,  1959,  Physics  and  Geology, 
McGraw  Hill,  NY.  - - 
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be  incorporated  into  the  calculations  of  this  report  (expressed  for  zero 
optical  depth):  2.8  x  10lc  (*a2)  s~\  2.8  x  10“  (ira2)  s'1,  2.8  x  10“  Ua2) 
s' 1 ,  and  0. 

Since  mesospheric  particles  may  possess  a  negative  or  positive  charge,  the 
question  arises  as  to  whether  an  electrostatic  potential  on  a  dust  particle 
would  modify  the  photodetachment  yield  from  that  dust.  Manka27  explored  the 
result  of  having  an  electrostatic  charge  on  the  surface  of  dust  and  the 
resultant  Coulomb  retarding  of  the  photoelectron  current.  He  was  specifically 
interested  in  dust  on  the  lunar  surface  and  showed,  for  example,  that  the 
current  density  of  photoelectrons  suffered  a  reauction  of  26  percent  in 
emitted  electrons  with  a  single  positive  charge.  fiur  calculations  are  not 
adjusted  for  this  small  erfect.  In  the  case  of  a  negative  potential,  the 
magnitude  of  the  photoolectron  current  was  not  modified.  That  is,  negatively 
charged  dust  will  provide  the  same  photoelectron  current  as  neutral  dust  and 
obey  the  same  wavelength  dependence.  Photoemissior.  wil’  not  be  enhanced. 

A  final  comment  is  necessary  about  the  differences  in  photoemission  between 
single  molecules  and  a  solid  composed  of  many  molecules:  It  was  considered 
possible  that  dust  could  photoemit  in  the  visible,  a  fact  suggested  by  analogy 
in  the  work  of  Woo  et  al  28  who  described  photoelectron  detachment  from 

the  0~  molecule  by  visible  wavelengths,  ’.his  idea  was  supported  by  Chanin  et 

al2“  who  showed  that  electrons  attaching  to  molecular  oxygen  to  form  0^  could 

possibly  enter  a  v  -  1  vibrational  state,  described  by  a  repulsive  potential 
energy  curve,  which  could  then  lead  to  autodetachment.  However,  in  the  case 
of  a  solid  dust  particle,  such  surface  metastable  states  of  high  energy  would 
not  exist  since  the  energy  of  the  impacting  electron  would  quickly  dissipate 
as  momentum  and  heat  through  neighboring  atoms.  The  "extra"  electron  causing 
the  unbalanced  negative  charge  on  a  dust  particle  would  be  indistinguishable 
from  any  other  electron  occupying  normal  bound  conduction  states. 
Photoemission  would  always  operate  against  the  work  function  of  the  material 
and  exhibit  the  same  cutoff  frequency  for  irradiation  as  that  for  an  uncharged 
surface.  Therefore,  it  would  not  be  likely  for  dust  to  photoemit  in  the 
vi  sibl e. 

4.3  Charge  Loading  on  Dust 

The  first  result  of  performing  the  calculation  of  the  effective  attachment 
coefficients,  and  K ^ ,  is  the  determination  of  the  fractions  of  the  dust 


27R.  H.  Manka,  "Plasma  and  Potential  at  the  Lunar  Surface,"  in  Photon  and 
Particle  Interactions  with  Surfaces  in  Space,  Proceedings  of  the' 6th  Esla~6 
"(yymposTumr  Noordwi The  ITetherl ands,  2^29  September  1972,  D.  Reidel 
Publishing  Co.,  Dordrecht-Hol 1  and/Boston  -  USA. 

2<S.  B.  Woo,  L.  M.  Branscomb,  and  E.  C.  Beaty.  1969,  "Sunlight  Photodetachment 
of  Ground  State  O^,"  J  Geophys  Res,  74:2933. 

2,L.  M.  Chanin,  A.  V.  Phelps,  and  M.  A.  Biondi,  1962,  "Measurements  of 
Attachment  of  Low  Energy  Electrons  to  Oxygen  Molecules,"  Phys  Rev,  128:219. 
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population  which  carry  a  particular  charge.  Table  2  lists  the  percentages  of 
the  total  dust  population  which  carry  given  numbers  of  electrons  for  dust  of 
radii  O.OOlym,  O.Olum,  and  0.1  um  under  various  ambient  electron  densities. 

The  tabulation  reveals  some  of  the  most  graphic  aspects  of  dust  behavior  in 
the  D  region:  For  example,  the  charge  for  dust  of  radius  O.Olum,  or  smaller, 
is  one  electron  per  dust  particle.  (These  results  are  similar  to  those  of 

figures  3  and  4  in  Parthasarathy .  ’ )  This  loading  holds  for  a  wide  '•anq<  of 
ambient  electron  number  densities  of  100  to  10,000  centimeter'’  an: 
photodetachment  values  of  0  to  10”  1  s"  !.  Under  simi  lar  ambient  cood’t' \ 
dust  sized  midway  between  O.Olum  to  0. m  heirs  a  more  or  less  stahV  ;.o 
of  2  to  3  electrons  which  shifts  to  '  ebstron  or  no  electrons  it  r  ' 
photodetachment  rates.  Oust  of  radius  m  presents  a  less  •. 1  *'■ a - c  * 

picture.  The  normal  loading  is  a  broad  di  si  •"l  :\>ti  on  of  5,  6.  or  ’  e  . 

which  is  displaced  to  positive  values  of  chargf'  where  photodetachment  . 

>1  s” 1  and  electron  densities  oelow  1000  cm'  ire  used.  In  each  of 
cases,  photodetachment  does  not  exert  an  in* ’oerce  until  it  exceeds  jbou* 
s'1. 

This  insensitivity  to  photodetachment  rates  below  a  certain  threshold  is 
important  result  since  a  vigorous  number-  for  dust  photodetachment  rates  hu" 
been  hard  to  obtain,  and  this  difficulty  had  boon  an  area  of  uncertainty  • 
chemical  models.  A  further  result  is  that  the  smaller  dust  sizes  will  <-xv  • 
in  the  D  -egion  with  only  a  single  negative  charge  if  they  are  charged  at  all. 

As  an  alternate  way  to  illustrate  the  competing  elements  of  this  situation, 
figure  3  shows  a  plot  of  photodetachment  freryency  for  moon  dust  and  ice  along 
with  electron  attachment  frequency  for  dust  rouged  Z  -  l,  0,  +1  and  frequency 
for  the  recombination  of  electrons  with  V"i‘  vd  hydrated  water  cluster  ions. 
An  ambient  electron  and  ion  density  o(  ’.  1  '  :  cm” '  has  been  used  for  tt-  c 
illustration. 

From  an  examination  of  figure  3,  a  number  of  points  are  apparent: 

a.  Positively  charged  dust  of  any  size  can  competitively  acquire 
electrons  when  compared  to  normal  electron-ion  recombination  processes. 
Uncharged  dust  of  radius  >  0.002uii  can  compete  effectively  for  electrons  even 
when  the  positive  ions  are  water  clusters  (as  is  the  case  for  altitudes  below 
80  km).  Negatively  charged  dust  can  only  compete  effectively  for  electrons  if 
the  particles  have  a  radius  >  0.02um.  Dust  with  radii  on  the  order  of  O.lutn 
could  attach  several  electrons,  but  because  of  the  nominally  low  number 
densities  of  O.ltm  radius  dust  ( v*  10  cm'3),  the  overall  electron  attachment  to 
dust  of  this  size  is  probably  of  secondary  importance  when  compared  to  the 
normal  electron-ion  recombination  processes,  a  point  already  noted  by 
Parthasarathy.  7  However,  if  dust  particles  of  much  smaller  sizes  were  present 
in  large  numbers,  as  is  most  likely  the  case,  then  dust  could  appreciably 
alter  the  ambient  electron  and  ion  equilibrium  densities.  The  precise  effect 


’R.  Parthasarathy,  1976,  "Mesopause  Dust  as  a  Sink  for  Ionization,"  J_Geophys 
Res,  81:2392. 
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TABLE  2.  PERCENT  OF  OUST  POPULATION  WITH  A  PARTICULAR 
NUMBER  OF  ELECTRONS  FOR  T  =  150°K* 


Dust  Particle  Radius  (urn) 


Ambient 


otodetachment 
Rate  (s'1) 


"Tnuur 


Percentage 


I'd,  DUST-ELECTRON  ATTACHMENT  FREQUENCY  ( 


EFFECTIVE  ATTACHMENT  COEFFICIENT,  K  „  &  K  *  El  FECTIVE  ATTACHMENT  COEFFICIENT 
FOR  DUST  OF  RADIUS  .001  um  FOR  DUST  OF  RADIUS  .01  um 


.“*•>)  SINaiOKUMOO  INaWHOVllV 
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effective  attachment  coefficients  of  O.OOlum  and 
to  electrons,  K®,  and  positive  ions,  KQ,  for 
electrons  and  ions  ranging  from  10  to  10“  cm-3. 


would  depend  on  the  rate  of  photoemi  ss’cn  of  the  electrons  from  the  neqativ*'1;. 
charged  dust. 


b.  During  twilight  or  nighttime,  smaller  particulates  'radius  ■. 
become  negatively  charged.  Hence,  in  largo  numbers,  small  dust  particles 
could  compete  effectively  with  io« -electron  ■r..  -'mbinatfon  and  provide  i  meg'"' 
of  maintaining  (rather  than  neutral i zing'  »ieg**ivn  charge.  Retorting  to  *r< 

right  side  of  equation  (b),  the  term  w'1'  nominate  the  tendency  tr  leal  ■' 

particles  negatively;  the  photnem-  •.-.  •  c  ,  *■  1  t‘”id  to 

may,  if  sufficiently  largo,  keep  *h  -.rticles  neut*- 3'  •  ve*<  • 

them  positive.  At  night,  or  fo>-  »•.->••; •  '  vr-.y  i-'v-i  < . 

will  dominate  and  drive  a  large  •  ■  •  ♦.’•••'  •  >--.i  c  1 '  3*0-  ,,  _ 

c.  Dust  which  is  ice  coated  photoemi ts  an  .cae-  of  magn’  tude  l>-<-  •  •' 

than  uncoated  dust.  This  difference  wed  1  or- an  that.  i..riuq  p.dar  ■ 
conditions  with  nocti lucent  clouds,  a  1.j>  ■■  number  of  electrons  co-i'  . 

attached  to  the  ice  coated  dust  than  at  otK  r  *;m;  s  of  the  v *-•>»-,  resul  f  in  3 

a  sharper  reduction  in  the  free  electron  >  tv 

4.4  Effective  Dust  Attachment  Coefficients 

Figure  4  shows  one  instructive  example  of  the  derived  tp  and  Kn  values  fc- 

O.OOlum  and  O.Olun  size  dust  plotted  against  electron  density  (Np)  anc  t.h.r  • 

photodetachment  values.  The  abscissa  for  figu*T  a  is  kj  (0r  M+  since,  by 

requirement  of  charge  neutrality,  they  must  be  held  equal  in 

calculations).  Immediately  apparent  from  th;s  graph  is  the  great  variabl'w 

of  Kp  and  Kp  with  changing  electron  and  positive  ion  densities  and 

photodetachment.  Specifically;  (a)  exceed  Kp  since  electrons  move  30U 

times  as  fast  as  positive  ions  and  hence  have  a  higher  collisional 

frequency.  (b)  Higher  values  of  photodetachment  rates.  vph,  result  in 

higher  Kp  values  since  photodetachment  terns  to  keen  a  particle  uncharged  (or 

positively  charged),  which  means  that  impacting  electrons  will  not  feel  a 

repulsive  force.  (c)  As  the  numbers  of  ambient  electrons  and  ions  increase. 

the  Kp  and  Kp  values  converge  to  a  common  asymptote  because  the  extreme 

ambient  densities  negate  the  effects  of  photodetachment.  However,  the  fact 

that  Kp  and  Kp  are  so  highly  dependent  upon  ambient  electron  and  ion  densities 

requires  that  they  be  entered  iteratively  into  computer  codes,  a  less 

efficient  procedure  than  using  a  set  of  kp  and  kp  for  the  individual  charging 
reactions  for  dust. 
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5.  THE  EFFECTS  OF  DUST  ON  D  REGION  CHEMISTRY  AND  CHARGE  BALANCE  THROUGH  THE 
USE  OF  A  FULL  CHEMICAL  CODE 

5.1  Runs  at  HO  km  with  Variable  Oust  Concentration  and  Photodetachment 
Frequencies 

The  results  of  the  previous  sections  have  been  organized  as  piecewise 
calculations  designed  to  gain  insight  into  the  behavior  of  dust  in  the 
mesosphere.  In  this  final  section,  dust  is  entered  as  a  full  partner  into  an 
electron-ion  reaction  set  to  quantify  the  effects  on  the  charge  balance.  The 
treatment  of  this  topic  is  not  exhaustive,  but  it  is  included  in  this  report 
to  establish  the  qualitative  nature  of  the  effect  of  dust  in  the  D  region. 

The  chemistry  reaction  set  is  that  of  C'AIRCHE^,  !  0  and  the  dust  attachment 

coefficients  used  are  -  i.?  1-4)  cm-’  :  for  0.01  urn  dust  (table  1).  Dust 

is  entered  in  several  concert  rut' one:  1,  '.00 ,  -''X,  --nd  10,000  cm*  5 ,  and  the 

resulting  values  in  electron  density.  NO*,  and  the  sum  of  the  hydrated  water 

cluster  ions  are  listed  in  table  3,  examples  1  through  4.  In  these  examples 

the  photodetachment  value  is  held  at  the  value  derived  from  Willis  et  al,25 
2.8  x  10'*  (-a*)  s"  \  and  the  dust  concentration  is  varied  by  4  orders  of 

magnitude.  A  second  set  of  examples,  5  through  7,  tests  the  results  using 
photodetachment  rates  which  vary  from  0.01  to  100  times  the  above  value. 

Examples  1  through  4  show  that  as  the  dust  population  increases  from  0  to  1000 
cm"'  and  10,000  cm"  ! ,  the  election  density  decreases  from  a  high  of  453  to  191 
and  54  cm"!,  respectively.  xty -  four  percent  of  dust  remains  uncharged  at 
high  concentrations  (10,000  cm":),  -nd  23  percent  is  negatively  charged. 
Besides  the  effect  of  reducing  the  electron  density,  another  major  change  is 
the  marked  increase  in  water  cluster  ion  population.  This  increase  is  due  to 
the  removal  of  electrons  which  would  recombine  with  the  cluster  ions  directly 
(that  is,  due  to  the  lower  ion-dust  ^combination  rate).  Total  charge 
neutrality  is  maintained  in  the  chemistry  of  these  runs,  although  only  a  small 
subset  of  positive  ion  species  entered  in  DAIRCHFM  is  listed  in  this  table. 
(For  example,  nitric  oxide  clusters  and  other  minor  positive  ion  species  are 
omi tted.  ) 

Examples  5  through  7  show  that  the  variation  in  photodetachment  is  significant 
in  modifying  the  free  electron  density.  High  phctodetachment  rates  drive  the 
dust  more  positive  and  hence  release  electrons.  cor  an  "increase  in 
photodetachment  of  4  orders  of  magnitude,  the  free  electrons  increase  from  118 
to  749  cm"3.  For  the  most  probable  value  of  photodetachment  (Willis'  value  of 
2.8  x  105  ta2  cm2  s~‘)2\  the  electron  density  increases  more  than  60  percent 
over  its  value  with  no  photodetachment. 


30D.  W.  Hoock  and  M.  G.  Heaps,  1978,  DAIRCHEM:  A  Computer  Code  to  Model 
.Ionizajtlon-Deionlzation  Processes  and  Chemi  stry  Tn  tKe  Middle  Atmosphere. 
Users  Manual,  OS  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile 
Range,  NM. 

I5R.  F.  Willis,  M.  Anderegg,  B.  Feuerbacher,  and  B.  Fitton,  "Photoemission  and 
Secondary  Electron  Emission  from  Lunar  Surface  Material,"  in  Photon  and 
Particle  Interactions  with  Surfaces  in  Space,  Proceedings  of  the  6th  TsTaF 
Symposium,  Noordwijk,  The  Netherlands,  26-2b-September  1972.  R.  J.  L.  Grand, 
Ed.,  D.  Reidel  Publishing  Co.,  Dordrecht-Holl and/Boston  -  USA. 
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In  each  set  of  test  runs  the  N0+  population  is  comparatively  stable.  In 
addition  to  the  reactions  listed  below,  a  set  of  attachment  reactions  of 
positive  ions  to  neutral,  negatively  charged,  and  positively  charged  dust  was 
also  used. 


Q  +  e  -  D' 

(8) 

D+  +  e  -  D 

(9) 

D-  +  hv  ♦  D  +  e 

(10) 

D+hv-D++e 

(11) 

5.2  Comparison  of  the  Effect  of  Four  Dust  Profiles  (r  =  O.Olum  and 
r  =  0.001pm)  upon  D-region  Ion  Chemistry. 

As  a  more  graphic  test  of  the  possibility  of  dust  modifying  the  electron 
profiles  of  the  0  region,  four  dust  distributions  were  included  in  the 
DAIRCHEM  code,  and  computations  were  made  at  a  series  of  altitudes.  As  a 
basis  of  comparison  for  the  modeling,  a  standard  type  of  observed  electron 
density  profile  from  daytime,  x  -  60°  at  Wallops  Island,  VA,  in  December  is 
presented.  Secondly,  the  electron  density  profile  computed  by  DAIRCHEM  with 
no  dust  included  is  plotted.  The  other  curves  are  computed  electron  densities 
using  0.01pm  dust  at  constant  concentrations  of  100  and  1000  cm"’  and  0.001pm 
dust  at  concentrations  of  1000  and  5000  cm-3.  These  values  bracket  the 
estimates  generally  accepted  for  dust  concentrations  in  that  region.  The 
values  are  tabulted  in  table  4.  A  photodetachment  value  of  2.8  x  lOM^a2)  s'* 
was  used  for  dust  for  all  runs. 


TABLE  4.  SELECTED  DUST  PROFILES 


Altitude  A  B  C  D  E 

(km)  No  Dust  0.01pm  0.001pm 


100 

0 

100 

1000 

1000 

90 

0 

100 

1000 

1000 

5000 

85 

0 

100 

1000 

1000 

80 

0 

100 

1000 

1000 

5000 

75 

0 

100 

1000 

1000 

70 

0 

100 

1000 

1000 

5000 

60 

0 

100 

1000 

1000 

The  dust 

profiles 

used  as 

i  nput  to  the 

code  represent 

a  plausible  range  of 

values. 

However, 

since  the  distribution 

of  dust  in 

the 

upper  mesosphere  is  a 

very  open 

question 

,  these 

profiles  were 

felt  to  be 

adequate  in  a  preliminary 

demonstration  of 

the  DAIRCHEM  results 

to  dust. 

particularly  the  charge 

27 


balance.  A  "step"  profile  of  dust  (i.e.,  a  sharp  ledqe  of  dust  at  a  given 

altitude)  has  been  conjectured  in  the  literature.  This  effect  may  easily  be 
seen  by  the  comparison,  at  any  altitude,  of  the  dust  case  with  the  no  dust 
case. 

The  results  of  these  inputs  to  DAIRCHEM  are  plotted  in  figure  5.  Curve  1  is 
the  Wallops  Island  observation  of  electron  density  for  December  with  y  = 

60°.  This  profile  shows  the  abrupt  change  in  slope  at  82  km  that  is  commonly 
observed  in  many  electron  density  profiles  in  the  D  region  under  quiet 

conditions.  Curve  2  is  the  computed  electron  density  profile  using  DAIRCHEM 

with  no  dust.  Note  the  divergence  of  these  two  curves  at  87  to  76  km.  Curves 
3  and  4  are  0.01  urn  dust  at  concentrations  of  100  cm' 5  and  1000  cm-3.  A 
reduction  in  electron  concentration  is  apparent  at  85  to  78  km,  particularly 
for  the  1000  cm' 3  distribution.  Below  78  km,  the  reduction  in  electron 
density  exceeds  the  observed  curve  1,  possibly  suggesting  that  dust  may  be 
confined  in  a  relatively  narrow  ledge.  For  example,  the  best  fit  using  0.01pm 
dust  would  be  in  excess  of  1000  cm'3  from  85  to  77  km  and  dropping  to  about 
200  cm' 3  by  70  km.  Curves  5  and  6,  O.OOlum  dust  at  1000  and  5000  cm'3, 
bracket  the  observed  values  of  electron  density  down  to  77  km.  Optimum 
concentration  for  0.001pm  dust  would  be  3000  to  4000  cm'5.  Relow  77  km  the 
reducing  action  of  this  dust  is  clearly  excessive,  thus  suggesting  values  near 
200  cm' 3  at  70  km. 

Although  these  calculations  were  run  at  each  5  km  in  altitude  and  with  a 
limited  number  of  cases,  the  results  clearly  show  that  a  judicious  choice  of 
dust  concentration  can  allow  the  modeled  electron  density  to  reproduce  the 
heretofore  "anomalous  bite-outs"  in  electron  density  often  observed  around  80 
km.  The  controlling  rule  is:  if  the  particle  concentration  equals  or  exceeds 
the  model  electron  density  computed  by  vapor  phase  reactions  alone,  such  as  in 
the  case  below  85  km,  a  reduction  to  observed  values  is  apparent.  Most  runs 
of  DAIRCHEM  indicate  the  best  fit  is  with  a  variable  concentration  of  dust 
with  altitude. 

6.  SUMMARY  AND  CONCLUSIONS 

The  purpose  of  this  report  has  been  to  examine,  through  a  series  of  piecewise 
calculations  and  comparisons,  the  importance  of  electron-dust  attachment  to  D 
region  chemistry,  particularly  in  competition  with  electron-ion 

recombination.  In  the  course  of  the  argument,  we  have  examined  the  preferred 
charge  loadings  on  a  population  of  dust  and  generated  two  different  (but 

related)  sets  of  coefficients  for  electron  or  positive  ion  dust  attachment 

6  + 

calculations.  The  discrete  attachment  coefficients,  k^  and  k^,  should 

generally  be  used  in  chemical  computer  codes  since  they  result  in  knowledge  of 
the  precise  charge  loading  of  individual  dust  particles. 

Dust  of  radius  less  than  O.Olim  acquires  only  one  electron  under  mesospheric 
conditions.  Oust  of  radius  0.1  urn  may  acquire  as  many  as  five  or  six  resident 
electrons,  and  in  this  size  range  photodetachment  does  not  perturb  the 
distribution  of  charge  loading  too  greatly.  The  size  of  dust  therefore  seems 
to  predetermine  a  unique  charge  loading  for  each  dust  particle  over  a  wide 
range  of  photodetachment  values.  An  accurate  measurement  of  dust  size  and 
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concentration  could  permit  one  to  reasonably  estimate  its  effect  upon  the 
chemistry  due  to  this  unique  property.  The  full  understanding  of  the  role  of 
dust,  however,  requires  the  use  of  the  ion  kinetic  reaction  scheme. 

This  report  ameliorates  some  uncertainty  in  photodetachment  rates  by  running  a 
set  of  calculations  of  DAIRCHEM  using  orders  of  magnitude  lower  and  higher 

than  the  value  used  by  previous  researchers  and  showing  resultant  stability  in 
the  computed  electron  densities.  Nevertheless,  the  previous  value,  obtained 
from  moon  dust,  seems  to  be  the  most  logical  for  use  in  further  D  region 
chemistry  runs.  Based  upon  the  limited  computations  made  for  this  study,  dust 
indeed  does  appear  to  compete  strongly  for  electrons,  particularly  at  dust 
concentrations  of  around  1000  cm-1. 

Figure  5  represents  the  most  advanced  set  of  calculations  of  this  report  and 
demonstrates  that  it  is  possible,  with  a  proper  choice  of  a  dust  profile  and 
size,  to  approximate  several  heretofore  "anomalous  features"  in  electron 
density  profiles. 

The  concentrations  of  dust  included  here  for  use  at  80  km  were  consistent  with 
reports  from  the  literature.  The  next  step  in  pursuing  this  topic  must  be  the 
analytic  derivation  of  a  dust  size  concentration  profile  starting,  most 
logically,  with  meteor  flux  estimates  at  150  km.  This  present  model  has  been 
run  strictly  for  daytime  chemistry  and  has  purposely  neglected  nighttime 
processes  and  the  attachment  of  negative  ions  to  dust  at  night.  Until  this 
other  important  step  is  taken  and  an  extensive  series  of  runs  is  made,  a 
description  of  the  role  of  dust  in  D  region  chemistry  will  remain  incomplete. 
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Department  of  Marine,  Earth 
and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27650 

Commander 

US  Army  Research  Office 

ATTN:  DRXRO-GS  (Or.  Leo  Alpert) 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Army  Research  Office 
ATTN:  DRXRO-PP  (Brenda  Mann) 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RD  (Dr.  K.  F.  Sterrett) 
Hanover,  NH  03755 

Commander/Di  rector 
US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRPEL-RG  (Mr.  George  Aitken) 

Hanover,  NH  03755 

Commander 

US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RG  (Mr.  Roger  H.  Berger) 
Hanover,  NH  C3755 

Commander 

US  Am\y  Armament  Research 
A  Development  Command 
ATTN:  DRDAR-AC  (Mr.  James  Greenfield) 

Dover,  NJ  07801 
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Commander 

US  Army  Armament  Research 
8  Oevc  lop.nert  Comma'' d 
ATTN:  DRPW-’SS  Shiny  »5?) 

Dover,  N.’  O’FOl 

Commander 

US  Army  Armament  Research 
&  Development  C-mmard 
ATTN:  DPCPM  "AW -••Cl  ?Mr.  Peteris  J arsons 
Dover,  S' 

Command0*' 

US  Army  Armament  Research 
8  Oeve1oemert  Command 
ATTN:  D :v.-,  u,  wa’dron' 
Dover.  N ' 

Deputy  Joint  ”-".:ject  v?."ager 
for  Navy/ J SMC  SAL  G7 

ATTN:  ORCO!>'-r,AWS-NV  'CCT  Joseph  MicelO 
Dover,  NJ  :)7gcl 

Commander/Di  rector 

US  Army  Combat  Surveillance  8  Target 
Acou’sit'or  Laboratory 
ATTN:  DELLS  - !  fMr.  David  tonginotti 1 
Fort  Monmculn,  NJ  O'  7  T 

Commander  /Di  rector 

US  Army  Combat  Surveillance  8  Targ°t 
Acquisition  Laboratory 
ATTN:  DCLCS-r:  '’■'r.  Sen  A.  Di  Campl  i } 
Fort  Monmouth,  NT  0 7 ? 0 2 

Commander 'D ; r°c cor 

■J S  Army  Combat  'y. r veil!  j nee  S  Target 
Ac qu i  f' t i on  L a :  o rate  ry 
ATTN:  DFLCS-P-S  'Mr.  Donald  L.  roian«) 
Fort  Monmnu*h  0’"03 

Director 

US  Army  Electronics  Technology  8 
Devices  Laboratory 
ATTN:  DELET-DD  (S.  Danko) 

Fort  Monmouth,  NJ  07703 

Project  Manager 
F1REFIN0ER/REMBASS 

ATTN:  DRCPM-FFR-TM  (Mr.  John  M.  BiaVo) 
Fort  Monmouth,  NJ  07703 


l  T. 

i:  '  c  ■ —ry  r  7  ect  r. >  i  r  n  Kes°arc° 

•  •  ■  oneet  L.>-ama  • 

A  '  o  L  -  LA  (Dr.  Walter  S.  H.  Afee) 
F-yrr  u  r.<<»,ah  N-l  077V- 


*  ‘  0"  Cir.  v-  C  '  y r  •'  %  Ac  *  ;  - 


v  ’?AvCC  Systems  Analysis  Activity 
V"*<:  AT  At- 1  OP  'Mr.  Louie  Dominguez) 
*,r-'”  Sands  Miss'le  Range,  NM  88003 


■  \  r*  t  »  •.  *  r\  *• 

yc  '  qa.0O';:  System?  Ara'ysis  Activity 
:Tii  -  '")■}  v  r,  william  J.  Leach) 

'■  '*  •  yarns  n:  -  Range,  NM  8300? 

?■'  •  ■:  .  co" 

-  'R.-V.-hC  Systems  Analysis  Activity 
ATAA-T  1?  '  i<!r.  Roger  F,  Willis) 
Tange,  NM  8800? 

"I recto^ 

-  'ice  f  Mlss'i:?  Electronic  Warfare 
ACTS:  ' '...L* -v-c '0  ‘’Dr.  Steven  Kovel) 

whs te  Ec-.r-d;  '*•  s  . Range.  NM  3800? 

Office  of  the  Test  Director 
Joint  Services  ED  GW  CM  Test  Program 
AtTN:  DRYDF-TD  (Mr.  Wei uOn  Findley! 
White  Sands  Missile  Range,  NM  88002 

Commander 

US  Army  White  Sands  Missile  Range 
ATTN:  STEWS-PT-AL  (laurel  B.  Saunders) 
Whit°  Sands  Missile  Range,  NM  88002 


40 


Commander 

US  Army  RSO  Coordinator 
US  Embassy  -  Bonn 
Box  165 

APO  New  York  09080 

Grumman  Aerospace  Corporation 
Research  Department  -  MS  A08-35 
ATTN:  John  E.  A.  Selby 
Bethpage,  NY  11714 

Rome  Air  Development  Center 
ATTN:  Documents  Library 
TSLD  (Bette  Smith) 

Griff iss  AFB ,  NY  13441 

Dr.  Roberto  Vagi io-Laurin 
Faculty  of  Arts  and  Science 
Dept,  of  Applied  Science 
26-36  Stuyvesant  Street 
New  York,  NY  10003 

Air  :'or-.;e  W-Ught  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN:  AF'WAL/AARI-3  (Mr.  Harold  Geltmacher) 
Wri  ght-Fatterson  AFB,  OH  45433 

Air  Force  Wright  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN:  AFWAL/AARI-3  (CPT  William  C.  Smith) 
Wri ght-Patterson  A^B,  OH  45433 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R  (CPT  James  M.  Watson) 
fort  Sill,  OK  73503 

Commandant 

US  " rmy  Field  Artillery  School 
ATTN:  ATSF-CD-MS 
Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R 
Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  NOAA  liaison  Officer 

(COR  Jeffrey  G.  Carlen) 

Fort  Sill,  OK  73503 


Commandant 

US  Army  Field  Artillery  School 
Morris  Swett  Library 
ATTN:  Reference  Librarian 
Fort  Sill,  OK  73503 

Commander 

Naval  Air  Development  Center 
ATTN:  '  -.e  301  (Mr.  George  F.  Eck) 

Warm! ns te-,  PA  18974 

The  University  of  Texas  at  El  Paso 
Electrical  Engineering  Department 
ATTN:  Dr.  Joseph  H.  Plerluissl 
E'  Paso,  TX  79963 

Corrmandant 

US  Army  Ai r  Defense  School 

ATTN:  ATSA-CD-SC-A  (CPT  Charles  T.  Thorn) 

Fort  31 iss,  TX  79916 

Commander 

HO.  7RAD0C  Combined  Arms  Test  Activity 
ATTN :  ATCAT-OP-Q  (CPT  Henry  C.  Cobb,  Jr.) 

Fort  Hoed,  TX  76544 

Commander 

HQ,  TRAQOC  Combined  Arms  Test  Activity 
ATTN :  ATCAT-SCI  (Dr.  Darrell  W.  Collier) 

Fort  Hood,  TX  76544 

Commander 

US  Army  Duoway  Proving  Ground 
ATTN :  STEDP-MT-DA-L 
Dug-way,  UT  84022 

Commander 

US  Armv  Ouaway  Proving  Ground 

ATTN:  S’ EDP-MT -DA-M  (Mr.  Paul  E.  Carlson) 

Dugway,  UT  84022 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-T  (Mr.  John  Trethewey) 

Dugway,  UT  84022 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-T  (Mr.  William  Peterson) 

Dugway,  UT  84022 
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')ef >'n -,r.  documentation  Center 
ATT'J  ;  0QC-TCA 
Cameron  Station  Dldg  5 
Alexandria,  V  A  22  314 
1? 


Ballistic  Miss'le  defense  Drogram  Office 
ATTN:  i)4CS-PMT  (Colonel  Harry  F.  Ennis) 
5001  Eisenhower  Avenue 
A 1  -?  x  a  n  d  r  i  a ,  V  A  2  2  3  3  3 

defense  Technica1  information  Center 
ArTN:  DDA-2  (Mr,  Carnes  c .  Chafer) 
Cameron  Station.  "Mq  5 
Alexandria.  VA  ? ? 3 1  -A 

Commander 

JS  Army  Materiel  Development 
&  Readiness  Command 

ATTN :  QRCRSi-EE  (Mr.  Albert  f.iambalvol 
5001  Eisenhower  Avenue 

Alexandria  V  A  2  ?  3  3  ; 

Commander 

'JS  Army  Materiel  Development 
3  Readiness  Command 
ATTN:  DRCi.OC  (Mr.  James  Render) 

5001  Eisenhower  Avenue 
Alexandria,  VA  223.33 

Oefensn  Advanced  Rsch  Projects  Agency 
ATTN:  Steve  Zakanvez 
1400  Wilson  (33  vd 
Arlington,  VA  22209 

Defense  Advanced  Rsch  Projects  Agency 
ATTN:  Dr.  James  Tegnelia 
1400  Wilson  B! vd 
Arlington,  VA  22209 

Institute  for  Defense  Analyses 
ATTN:  Mr.  Lucien  M.  Biberman 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Institute  for  Defense  Analyses 
ATTN:  Dr.  Ernest  Bauer 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Institute  for  Defense  Analyses 
ATTN:  Dr.  Hans  G.  Wolfhard 
400  Army-Navy  Drive 
Arlington,  VA  22202 


System  Planning  Corporation 
ATTN:  Mr .  Daniel  rriedman 
1500  W' 1  son  Boulevard 
Arlinqton,  VA  22209 

System  Planning  Corporation 
ATtN:  COL  Hank  Shelton 
1500  Wilson  Boulevard 
Arlington,  VA  22209 

C S  a -my  Intelligence  ft  Security  Command 
A'r’’-.:  Edw'r  Speakman.  Scientific  Advisr- 
A-’inotor  Via "1 1  Station 
A-l  mgton,  VA  22212 

Commander 

Army  Operational  Test 
ft  :  val  t;at  ’  on  Agency 
ATTN:  CB7E-FD  (Mr.  Floyd  I.  Hill) 

5600  Columbia  Pike 
Falls  Church,  VA  22041 

Commander  and  Director 
'JS  Army  Engineer  Topographic  Laboratories 
ATTN:  ETL-GS-A  (Mr.  Thomas  Nei dringhaus 1 
Fort  Bel  voir,  VA  22060 

Director 

'JS  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN:  DELNV-L  (Dr.  Rudolf  G.  Buser) 

Fort  Bel  voir,  VA  22060 

Di rector 

US  Army  Night  Vision  ft 
Electro-Optics  Laboratory 
ATTN:  DELNV-L  (Dr.  Robert  S.  Rodhe) 

Fort  Bel  voir,  VA  22060 

Di rector 

US  Army  Night  Vision  ft 
Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Mr.  Joseph  R.  Moulton) 
Fort  Bel  voir,  VA  22060 

Director 

US  Army  Night  Vision  ft 
Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Luanne  P.  Obert) 

Fort  Bel  voir,  VA  22060 
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Di rector 

'J S  Army  Night  Vision 
A  Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Mr.  Thomas  W.  Cassidy) 
Fort  Bel  voir,  VA  22060 

Oi rector 

US  Army  Night  Vision  A 
Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Mr.  Richard  J.  8ergemann) 
Fort  Bel  voir,  VA  22060 

Oi rector 

US  Army  Night  Vision  ft 
Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Or.  James  A.  Patches) 

Fort  Del  voir,  VA  22060 

Commander 

US  Army  Training  ft  Ooctrine  Command 

ATTN:  ATCD-AN 

cort  Monroe,  VA  23651 

Commander 

US  Army  Traininq  ft  Doctrine  Command 

ATTN:  ATCD-AN-M 

Fort.  Monroe,  VA  23651 

Commander 

US  Army  Training  ft  Doctrine  Command 
ATTN:  ATCO-F-A  (Mr.  Chris  O'Connor,  Jr.) 
Fort  Monroe,  VA  23651 

Commander 

IJS  Army  Training  ft  Doctrine  Command 
ATTN:  ATCD-IE-R  (Mr.  David  M.  Inqram) 

Fort  Monroe,  VA  23651 

'.ommander 

US  Army  Training  ft  Doctrine  Command 
ATTN:  ATCD-M-I/ATCD-M-A 
::ort  Monroe,  VA  23651 


Department  of  the  Air  Force 
HQS  5  Weather  Wing  (MAC) 

ATTN:  5  WW/QN 

Langley  Air  rorce  Base,  VA  23655 
Commander 

US  Army  INSCOM/Quest  Research  Corporation 
ATTN:  Mr.  Donald  'Wilmot 
6845  Elm  Street,  Suite  407 
McLean,  VA  22101 

General  Research  Corporation 
ATTN:  Dr.  Ralph  Zirkind 
7655  Old  Springhouse  Road 
McLean,  VA  22102 

Science  Applications,  Inc. 

8400  Westpark  Drive 
ATTN:  Dr.  John  E.  Cockayne 
McLean,  VA  22102 

US  Army  Nuclear  ft  Chemical  Agency 
ATTN:  MONA-WE  (Dr.  John  A.  Berberet) 

7500  Backlick  Road,  Bldg  2073 
Springfield,  VA  22150 

Di rector 

US  Army  Siqnais  Warfare  Laboratory 
ATTN:  DELSW-EA  (Mr.  Douglas  Harkleroad) 
Vint  Hill  Farms  Station 
Warrenton,  VA  22186 

Di rector 

US  Army  Signals  Warfare  Laboratory 
ATTN:  DELSW-OS  (Dr.  Royal  H.  Rurkhardt) 
Vint  Hill  Farms  Station 
Warrenton,  VA  22186 

Commander 

US  Army  Cold  Regions  Test  Center 
ATTN:  3TECR-TD' (Mr,  Jerold  Barger) 

APO  Seattle,  WA  98733 


Commander 

US  Army  Training  A  Doctrine  Command 
ATTN:  ATDOC-TA  (Dr.  Marvin  P.  Pastel) 
Fort  Monroe,  VA  23651 

Department  of  the  Air  Force 
OL-I,  AWS 

Fort  Monroe,  VA  23651 


HQDA  (SAUS-OR/Hunter  M.  Woodall,  Jr./ 
Dr.  Herbert  K.  Fallin) 

Rm  2E  614,  Pentagon 
Washington,  DC  20301 

COL  Elbert  W.  Friday,  Jr. 

OUSORE 

Rm  30  129,  Pentagon 
Washington,  DC  20301 
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Defense  Communications  Agency 
'echnical  Library  Center 

Code  ?;>? 

Washington,  nr  20305 
D i rector 

defense  Nuclear  Agency 

ATTN:  Technical  Library  (Mrs.  Betty  Fox) 
Washington,  DO  20305 

Di rector 

be  Tense  Nuclear  Agency 
ATTN:  R.AAE  (Or.  Carl  Fitz) 

Wash ; ngton,  OC  20305 

Oi rector 

Defense  Nuclear  Aqency 

ATTN;  SPAS  (Mr.  Donald  J.  Kohler) 

Washington,  DC  20305 

Defense  Intelligence  Agency 
ATTN:  DT/AC  (HC  Robert  Poplawski) 
Washington,  OC  20301 

HQOA  (DAMA-ARZ-O/Or.  Verderame) 
Washington,  OC  20310 

HQOA  (DAMI-ISP/Mr.  Beck) 

Washington,  DC  20310 

Department  of  the  Army 
Deputy  Chief  of  Staff  for 
Operations  and  Plans 
ATTN:  DAMO-RQ 
Washington,  DC  20310 

Department  of  the  Army 
Director  of  Telecommunications  and 
Command  and  Control 
ATTN:  DAMO-TCZ 
Washington,  DC  20310 

Department  of  the  Army 

Assistant  Chief  of  Staff  for  Intelligence 

ATTN:  DAMI-TS 

Washington,  DC  20310 

HQOA  (DAEN-ROM/Qr.  de  Percin) 

Casimir  Pulaski  Buildinq 
20  Massachusetts  Avenue 
Room  6203 

Washington,  OC  20314 


National  Science  Foundation 
Division  of  Atmospheric  Sciences 
ATTN:  Dr.  Eugene  W.  Rierly 
1800  G.  Street,  N.W. 

Washington,  OC  20550 

Director 

Naval  Research  Laboratory 

ATTN:  Code  4320  (Dr.  Lothar  H.  Ruhnke) 

Washington,  DC  20375 

Command' no  Officer 

Nava1  Research  Laboratory 

A^TN:  Code  600°  (Dr.  John  vacCailjm,  ’r 

Washington,  DC  20375 

Commanding  Officer 

Naval  Research  Laboratory 

ATTN:  Code  6530  (Mr.  Raymond  A.  Patten'' 

Washington,  DC  20375 

Commanding  Officer 

Naval  Research  Laboratory 

ATTN:  Code  6533  (Dr.  James  A.  Dowling) 

Washington,  DC  20375 
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ATMOSPHERIC  SCIENCES  RESEARCH  REPORTS 

1.  Lindberg,  J.  D.  'An  Improvement  to  a  Method  for  Measuring  the  Absorption 

Coefficient  of  Atmospheric  Dust  and  other  Strongly  Absorbing 
Powders,"  ECOM-5565,  July  197S. 

2.  Avara,  Elton  P. ,  "Mesoscale  Wind  Shears  Derived  from  Thermal  Winds," 

EC0M-5566,  July  1975. 

3.  Gomez,  Richard  B. ,  and  Joseph  H.  Pierluissi,  "Incomplete  Gamma  Function 

Approximation  for  King's  Strong-Line  Transmittance  Model," 
ECOM-5567 ,  July  1975. 

4.  Blanco,  A.  J.t  and  B.  F.  Engebos,  "Ballistic  Wind  Weighting  Functions  for 

Tank  Projectiles,"  ECOM-5568,  August  1975. 

5.  Taylor,  Fredrick  J.,  Jack  Smith,  and  Thomas  H.  Pries,  "Crosswind 

Measurements  through  Pattern  Recognition  Techniques,"  ECOM-5569, 
July  1975. 

6.  Walters,  D.  L. ,  "Crosswind  Weighting  Functions  for  Direct-Fire 

Projectiles,"  EC0M-5570,  August  1975. 

7.  Duncan,  Louis  D. ,  "An  Improved  Algorithm  for  the  Iterated  Minimal 

Information  Solution  for  Remote  Sounding  of  Temperature,"  ECOM-5571, 
August  197b. 

8.  Robbiani,  Raymond  L. ,  "Tactical  Field  Demonstration  of  Mobile  Weather 

Radar  Set  AN/ TPS -4 1  at  Fort  Rucker,  Alabama,"  ECOM-5572,  August 
1975. 

9.  Miers,  B. ,  G.  Blackman,  D.  Langer,  and  N.  Lorimier,  "Analysis  of  SMS/GOES 

Film  Data,"  ECOM-5573,  September  1975. 

10.  Manquero,  Carlos,  Louis  Duncan,  and  Rufus  Bruce,  "An  Indication  from 

Satellite  Measurements  of  Atmospheric  CO2  Variability,"  ECOM-5574, 
September  1975. 

11.  Petracca,  Carmine,  and  James  D.  Lindberg,  "Installation  and  Operation  of 

an  Atmospheric  Particulate  Collector,"  EC0M-5575,  September  1975. 

12.  Avara,  Elton  P. ,  and  George  Alexander,  "Empirical  Investigation  of  Three 

Iterative  Methods  for  Inverting  the  Radiative  Transfer  Equation," 
EC0M-5576,  October  1975. 

13.  Alexander,  George  D. ,  "A  Digital  Data  Acquisition  Interface  for  the  SMS 

Direct  Readout  Ground  Station  -  Concept  and  Preliminary  Design," 
EC0M-5577,  October  1975. 

14.  Cantor,  Israel,  "Enhancement  of  Point  Source  Thermal  Radiation  Under 

Clouds  in  a  Nonattenuating  Medium,"  EC0M-5578,  October  1975. 
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15.  Norton,  Colburn,  and  Glenn  Hoidale,  "The  Diurnal  Variation  of  Mixing 

Height  by  Month  over  whi'e  Sands  Missile  Range,  NM,"  ECOM-5579, 
November  1975. 

16.  Avara,  Elton  P.  ,  "On  the  Spectrum  Analysis  of  binary  Data,"  EC0M-5580, 

November  1975. 


17.  Taylor,  Fredrick  J.,  Thomas  H.  Pries,  and  Chao-Huan  Huang,  "Optimal  Wind 
Velocity  Estimation,"  ECOM-5581,  December  1975. 


13.  Avara,  Elton  P. ,  "Some  Effects  of  Autocorrelated  and  Cross-Correlated 
Noise  on  the  Analysis  of  Variance,"  ECOM-5582,  December  1975. 


19. 


Gillespie,  Patti  S. ,  R.  L.  Armstrong,  and  Kenneth  0.  White,  "The  Spectral 
Character i sti cs  and  Atmospheric  COo  Absorption  of  the  Ho+^:YLF  laser 
at  2.05um,"  ECOM-5583,  December  1975. 


20.  Novlan,  David  J.,  "An  Empirical  Method  of  Forecasting  Thunderstorms  for 

the  White  Sands  Missile  Range,"  ECOM-5584,  February  1976. 

21.  Avara,  Elton  P.,  "Randomization  Effects  in  Hypothesis  Testing  with 

Autocorrelated  Noise,"  EC0M-5585,  February  1976. 

22.  Watkins,  Wendell  R.,  " Improvements  in  Long  Path  Absorption  Cell 

Measurement,"  EC0M-5586,  March  1976. 


23.  Thomas,  Joe,  George  D.  Alexander,  and  Marvin  Dubbin,  "SATTEL  -  An  Army 

Dedicated  Meteorological  Telemetry  System,"  ECOM-5587,  March  1976. 

24.  Kennedy,  Bruce  W.,  and  Delbert  Bynum,  "Arny  User  Test  Program  for  the 

RDTSE-XM-75  Meteorological  Rocket,"  ECOM-5588,  April  1976. 

25.  Barnett,  Kenneth  M. ,  "A  Description  of  the  Artillery  Meteorological 

Comparisons  at  White  Sands  Missile  Range,  October  1974  -  December 
1974  ('PASS’  -  Prototype  Artillery  [Meteorological]  Subsystem)," 
E COM-5589,  April  1976. 

26.  Miller,  Walter  B. ,  "Preliminary  Analysis  of  Fall-of-Shot  From  Project 

'PASS',"  EC0M-5590,  April  1976. 

27.  Avara,  Elton  P. ,  "Error  Analysis  of  Minimum  Information  and  Smith's 

Direct  Methods  for  Inverting  the  Radiative  Transfer  Equation," 
EC0M-5591 ,  April  1976. 

28.  Yee,  Young  P. ,  James  D.  Horn,  and  George  Alexander,  "Synoptic  Thermal 

Wind  Calculations  from  Radiosonde  Observations  Over  the  Southwestern 
United  States,"  ECOM-5592,  May  1976. 
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29.  Duncan,  Louis  D. ,  and  Mary  Ann  Seagraves,  "Applications  of  Empirical 

Corrections  to  NOAA-4  VTPR  Observations ECOM-5593,  May  1976. 

30.  Miers,  Bruce  T. ,  and  Steve  Weaver,  "Applications  of  Meteorological 

Satellite  Data  to  Weather  Sensitive  Army  Operations,"  ECOM-5594,  May 
1976. 

31.  Sharenow,  Moses,  "Redesign  and  Improvement  of  Balloon  ML-566,"  ECOM-5595, 

June  1976. 

32.  Hansen,  Frank  V. ,  "The  Depth  of  the  Surface  Boundary  Layer,"  ECOM-5596, 

June  1976. 

33.  Pinnick,  R.  G. ,  and  E.  B.  Stenmark,  "Response  Calculations  for  a 

Commerical  Light-Scatteri ng  Aerosol  Counter,"  ECQM-5597,  July  1976. 
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